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Abstract. Consider the unit triangular lattice in the plane with origin O , drawn so 
that one of the sets of lattice lines is vertical. Let l and V denote respectively the 
vertical and horizontal lines that intersect O. Suppose the plane contains a pair of 
triangular holes of side length two, distributed symmetrically with respect to l and l', 
and oriented so that both holes point toward O. Unit rhombus tilings of three different 
regions of the plane are considered, namely: tilings of the entire plane; tilings of the 
half plane that lies to the left of l (where l is considered a free boundary, so unit 
rhornbi are allowed to protrude halfway across it); and tilings of the half plane that 
lies just below the fixed boundary l'. Asymptotic expressions for the interactions of 
the triangular holes in these three different regions are obtained, providing further 
evidence for Ciucu’s ongoing program that seeks to draw parallels between gaps in 
dimer systems on the hexagonal lattice and certain electrostatic phenomena. 


1. Introduction 

The study of interactions between holes (or gaps) in dimer systems was established by 
Fisher and Stephenson [T5] in 1963 with their results concerning interactions of holes in 
dimer systems on the square lattice. Three specific types of interaction were presented: 
the interaction between a pair of dimer holes (that is, the interaction of two fixed dimers 
within a dimer system); the interaction between a pair of non-dimer holes (that is, the 
interaction between two fixed monomers in a dimer system); and the interaction of a 
dimer hole with a constrained boundary. 

Kenyon [20] generalised the first of these results to an arbitrary number of dimer 
holes in dimer systems on both the square and hexagonal lattices. In [22] Kenyon, 
Okounkov and Sheffield further extended this work to general planar bipartite lattices. 
The interaction of a monomer with a straight line boundary in the square lattice is 
also due to Kenyon [21] Section 7.5], while the interaction of a family of holes with a 
constrained straight line boundary on the hexagonal lattice is tackled by Ciucu in [5] 
Theorem 2.2], It is worth noting here that since the triangular lattice arises as the 
“dual” of the hexagonal lattice, dimer coverings of subgraphs of the hexagonal lattice 
are often much more easily thought of in terms of unit rhombus tilings of sub-regions 
of the triangular lattice in the plane (so monomers correspond to unit triangles and 
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dimers correspond to unit rhombi). A dimer covering of the entire hexagonal lattice 
then corresponds to a tiling of the entire plane by unit rhombi (sometimes referred to 
as a sea of unit rhombi). 

Ciucu has also studied various types of interactions between non-dimer holes in dimer 
systems on the hexagonal lattice (equivalently, interactions between triangular holes 
within a sea of unit rhombi [3] [3j [2] [6] [8j ), thereby establishing close analogies to cer¬ 
tain electrostatic phenomena. Ciucu conjectures that interactions between holes in 
dimer systems on the hexagonal lattice are governed by the laws of two dimensional 
electrostatics (that is, Coulomb’s laws). More specifically, the conjecture states that 
by taking the exponential of the negative of the electrostatic energy of the two dimen¬ 
sional system of physical charges obtained by considering each hole as a point charge of 
size and magnitude specified by a statistical q, one may recover (up to a multiplicative 
constant) asymptotic expressions for the dimer-mediated interactions of holes on the 
hexagonal lattice. This conjecture remains wide open, though it has been shown to 
hold for two fairly general families of holes. 

In light of this, further types of interaction have been studied. Recently in [13], Ciucu 
and Krattenthaler determined the interaction between a left-pointing triangular hole 
of side length two and a vertical free boundary that borders a sea of unit rhombi on 
the left half of the plane. It would appear that this result (which inspired the present 
paper) is the first treatment of such an interaction in the literature. Theorem !2.2l below 
is a direct analogy of the interaction presented in pj3] in the sense that here the hole has 
been “flipped” and instead points toward the free boundary. Curiously, Theorem 12.21 
agrees with that of [13] only up to a multiplicative constant. 

As well as Theorem 12.21 this paper establishes two other types of interaction between 
a pair of triangular holes of side length two that point towards each other in dimer 
systems (that is, a pair of triangles oriented so that one left-pointing triangle lies directly 
to the right of one right-pointing triangle). These interactions are: the interaction 
between such a pair of holes that lie within a sea of unit rhombi (Theorem 12.ip : and 
the interaction between such a pair of holes that lie on a fixed boundary of a sea of 
unit rhombi lTheorem l2.3p . All three of these results agree with the original conjecture 
of Ciucu and so serve as yet more evidence that interactions between holes in dimer 
systems are governed by the laws of two dimensional electrostatics. A brief outline of 
the paper follows. 

Section [2] discusses the correspondence between dimer systems on the hexagonal lat¬ 
tice and rhombus tilings of hexagons on the triangular lattice, and also states precisely 
the family of holey hexagons (that is, hexagons containing triangular holes) that are 
of particular interest. Certain symmetry classes of rhombus tilings of a holey hexagon 
correspond precisely to rhombus tilings of smaller sub-regions of the same hexagon. The 
three different types of interaction that are the focus of this work arise from considering 
tilings of these sub-regions as they become infinitely large. These sub-regions are de¬ 
fined in Section [2] along with the corresponding interactions (or correlation functions). 
Asymptotic expressions for these interactions are also stated here fTheorem l2.ll Theo¬ 
rem [272] and Theorem 12.3p . which together comprise the main results. Proposition 12.41 


1 For a triangular hole t on the unit triangular lattice, q(t) is defined to be the signed difference 
between the number of left and right-pointing unit triangles that comprise t. 
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a factorisation result that is frequently used throughout this paper- is also stated in 
this section. 

Section [H] briefly discusses the relationship between rhombus tilings of hexagons and 
plane partitions. This is followed by four exact enumeration formulas that count tilings 
of sub-regions of holey hexagons (Theorem 13.11 Theorem 13.21 Theorem 13.31 and Theo¬ 
rem 13.411 . Proofs of these theorems follow in Sections [4] and [51 

In Section [4] the well-known bijection between rhombus tilings and non-intersecting 
lattice paths is combined with existing results due to Stembridge [31] and Lindstrom- 
Gessel-Viennot [IS] in order to express the number of tilings of different regions as 
either a Pfaflian (Proposition [O]) or a determinant (Proposition 14.41 Proposition 14.51) . 
In Section ED the matrices defined in Proposition 14.31 are manipulated using standard row 
and column operations in such a way that the Pfaflian of each matrix may be expressed 
as a determinant of a much smaller matrix. The unique Lf/-decompositions of these 
smaller matrices, along with those from Proposition 14.41 and Proposition 14.51 are then 
stated and proved (Theorem 15.21 Lemma 15.31 Theorem 15.41 and Theorem 15.51) . from 
which the exact enumeration formulas in Section [3] immediately follow. 

Section [6] contains proofs of the main results. These follow from inserting the enu¬ 
meration formulas from Section [3] into the correlation functions defined in Section [21 
Asymptotic expressions for the three types of interaction are then derived by studying 
these correlation functions as the distance between the holes (or the distance between 
the hole and the free boundary) becomes very large. 


2. Set-Up and Results 


Consider the hexagonal lattice in terms of its dual 3F (that is, the planar triangular 
lattice consisting of unit equilateral triangles drawn so that one of the sets of lattice 
lines is vertical). Then monomers on correspond to unit triangles on 2F and dimers 
on correspond to unit rhombi on SF (a dimer on Ftf? is equivalent to joining two 
unit triangles that share exactly one edge on £F). It follows that dimer coverings of 
regions of that contain gaps correspond exactly to rhombus tilings of regions of 2F 
that contain triangular holes. 

Fix some origin O in SF and denote by R a ,b the hexagonal region centred at O with 
side lengths a, b, a, a, b, a (going clockwise from the southwest edge). Suppose T is 
a set of fixed unit triangles contained within the interior of H ab , where the distances 
between the triangles are parametrised by k. Denote by H a \ T the hexagonal region 
H a ^ with the set of triangles T removed from its interior (such regions will often be 
referred to as holey hexagons, a term coined by Propp in [28]). Suppose further that 
£ > 0 is real such that b ~ £a is integral. Then the interaction (otherwise known as the 
correlation function, or simply the correlation ) between the holes in T is defined to be 


w(fc;£) 


lim 

a—> oo 


M{H a)b ) 


where M(R) denotes the number of rhombus tilings of the region R. Different types of 
interactions may be obtained by replacing H ab \ T and H a b with specific sub-regions of 
H a ,b \ T (and H a b respectively) in the above formula. Note that in order for a rhombus 
tiling of H ay b to exist (and hence for the above definition to make sense), it must be the 
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Figure 1. Left: the holey hexagon -H 7 , 5 \(> 4 U< 4 ). Right: a small section 
of the plane of unit triangles given by H ab \{> 9 U<g) as a tends to infinity. 


case that ^2 t€T q(t) = 0, where q(t) is the statistic defined in Section [TJ The triangular 
holes and sub-regions of particular interest are defined below. 

For some non-negative integer k, let >k denote the right-pointing equilateral triangle 
of side length two at lattice distance k to the left of the origin such that the centre 
of its vertical side lies on the horizontal symmetry axis of H a Define <k analogously, 
that is, denotes the left-pointing triangle of side length two with its vertical edge at 
lattice distance k to the right of the origin. Then H ab \ (>*. LM*,) denotes the hexagonal 
region H a jb from which two inward pointing triangular holes of side length two have been 
removed that are symmetrically distributed with respect to the vertical and horizontal 
symmetry axes of H a>b . Figure Q] (left) shows such a region for a = 7, b = 5 and k = 4. 

Remark 2.1. Under this construction the parity of a and b determine the parity of k, 
that is, if a and b are of the same parity then k is even, otherwise k is odd. 


Consider H ab \ (>k U <&) as a is sent to infinity. It should be clear that this gives the 
entire plane of unit triangles containing two triangular holes with their vertical sides 
at lattice distance 2k apart, a portion of which may be seen in Figure [Tj (right). Hence 
the correlation function 


u H {k-£) 


lim 

a—>oo 


M(H a ,b \ (>k U <k )) 
M(H aib ) 


is the interaction (in terms of k) between a pair of triangular holes of side length 2 
oriented so that they point directly at each other within a sea of unit rhombi. 


Theorem 2.1. As k becomes very large the interaction defined above, ujuik-,^), is 
asymptotically 


UH{k;£) ~ 




Denote by V a ^ \ >k the left half of H a b \ (> k U <*,) whose rightmost boundary is the 
vertical line l that intersects the origin. If l is a free boundary (that is, unit rhombi are 
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FIGURE 2 . Left: a vertically symmetric rhombus tiling of Hq^ \ (>4U<4). 
Centre: the corresponding tiling of Vq^ \ > 4 . Right: a small section of 
the left half of the plane of unit triangles given by V a ^ \ > 2 o as a tends to 
infinity. 


permitted to protrude across /), then the number of rhombus tilings of V a}b \>k is clearly 
equal to the number of vertically symmetric tilings of H ab \ (> k U <k). Figure [2] (left) 
shows such a tiling for H e 6 \ (> 4 U < 4 ) along with the corresponding tiling of Ve )6 \ > 4 
(centre). 

As a tends to infinity, the region V a ^ \ gives the left half plane of unit triangles, 
constrained on the right by /, and containing a right-pointing triangular hole with 
its vertical side at lattice distance k to the left of l (see Figure [21 right). Hence the 
correlation function 


cny(fc;0 


lim 

a—>• 00 


M(V a ,b \ >fc) 
M(V a , b ) 


gives the interaction (in terms of k ) between a right-pointing triangle of side length two 
and a vertical free boundary that borders a sea of unit rhombi that tile the left half 
plane. 


Theorem 2.2. As k becomes very large the interaction defined above, 
asymptotically 


uv{k]Q ~ 


yW±A 

2nk 


(£ + 1 ) 

2 


2(fc—1) 


Uv{k)€), 


is 


Consider now the two regions obtained by cutting H a b \ U <*,) along the zig-zag 
line that lies just below its horizontal symmetry axis. Denote the upper and lower 
portions of H atb \ (> k U < fc ) with Hfi b \ (t> fc U < fc ) and H~ b \ (> fe U < fe ) respectively. Figure El 
(left) shows the upper and lower regions of H 7 8 \ (t > 5 U < 5 ). If b is even then rhombus 
tilings of H~ b \(> k U<fc) correspond to horizontally symmetric tilings of the entire region 
Ha,b \ (>fc U <h). 

Now consider H~ b \ (> k U < k ) as a is sent to infinity. Clearly one obtains the lower 
half plane of unit triangles, bounded above by a fixed boundary on which lie two 
triangular holes with vertical sides at lattice distance 2k apart- see Figure [3] (right). 
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FIGURE 3. Left: The regions Hf 8 \ (> 4 U < 4 ) (upper) and Hf 8 \ (> 4 U < 4 ) 
(lower). Right: a small section of the lower half of the plane of unit 
triangles given by H~ b \ (>7 U < 7 ) as a tends to infinity. 


The correlation function 


wjj-(fc;0 


lim 

a—>00 




then gives the interaction (in terms of k) of a pair of triangular holes of side length 
two that point directly towards each other and are positioned along the fixed upper 
boundary that borders unit rhombus tilings of the lower half plane. 


Theorem 2.3. As k becomes very large the correlation function defined above, u)H-{k ; £), 
is asymptotically equivalent to that ofujy(k]f) from Theorem \2.°A 


Observe that for £ 7 ^ 1, Theorem 12.11 and Theorem 12.21 give distorted dimer statis¬ 
tics. Indeed, these expressions blow up or vanish exponentially for £ > 1 and £ < 1 
respectively. So when considering the three types of interaction detailed previously, it 
makes sense to consider only regular hexagons, that is, hexagons with all sides the same 
length. I 11 light of this the above asymptotic expressions may be re-packaged in the 
following way: 


3 

u v(k]l) ~ . 

47rd(>fc, O) 

UH(k;1} ~ 


( 2 . 1 ) 

( 2 . 2 ) 


where d(a, b ) denotes the Euclidean distance between a and b. 

The first result above is analogous to that of Ciucu and Krattenthaler [I3j, which 
presents an asymptotic result for the correlation of a triangular hole that has been 
“flipped” and instead points away from the vertical free boundary l. Within an infinitely 
large sea of rhombi one could be forgiven for thinking that the orientation of such a 
hole would have no effect on the interaction between the hole and the free boundary, 










































































THREE INTERACTIONS OF HOLES IN TWO DIMENSIONAL DIMER SYSTEMS 


7 


yet this result shows that such an interaction differs by a factor of 3. Any intuitive 
reason as to why this phenomenon should occur remains a complete mystery. 

Both of these results are entirely in keeping with the original conjecture of Ciucu |4j, 
Conjecture 1]. The asymptotic expression for ojuik^ff) may be obtained (up to a mul¬ 
tiplicative constant) by considering the exponential of the negative of the electrostatic 
energy of the system of physical charges obtained by viewing the triangular holes as 
point charges, where one has signed magnitude +2 and the other —2 according to the 
definition of the statistic q in Section [0 In a similar way the asymptotic expression 
for (jj v {k ;£) shows that a triangular gap is attracted to a vertical free boundary l in 
precisely the same way that an electrical charge is attracted to a straight line conductor 
when placed near it. Hence (12.11) and (12. 2\) above further support the hypothesis that 
interactions of gaps in the hexagonal lattice are governed by Coulomb’s laws for two 
dimensional electrostatics. 

This section concludes with an important proposition concerning the enumeration of 
tilings of the entire region H ab \ (> k U <*,). It follows from the Matchings Factorization 
Theorem [7] by employing arguments that are almost identical to those of Ciucu and 
Krattenthaler [12j . 


Proposition 2.4. The total number of tilings of the region H ab \ (>*. U < k ) has the 
factorisation 

M(H a , b \ (> fe U <*)) = M(Hf b \ ( > k U <*)) • M w {H+ b \ (> k U <*)), 

where M(H) denotes the number of tilings of the region H and M w (H+ b \ (>*. U <&)) 
denotes the weighted count of the region H+ h \ (>*. U <*,) where the unit rhombi that lie 
along its bottom edge carry a weight of 2. 


Remark 2.2. From now on, denote by Hf b \ (>*. U < k ) the region Hff b \ (> k U <*,) with 
the added condition that the unit rhombi lying along its bottom edge have a weight of 
2, so that M w (H+ b \ (> fc U < k ) = M(H+ b \ (> k U < k )). 


This proposition implies that rhombus tilings of H a b \ U <^) may be enumerated 
by counting separately the number of tilings of Hf b \ U <*,) and H~ b \ (> k U < k ) and 
taking their product. If b is an even integer then tilings of the latter region correspond 
to horizontally symmetric tilings of the entire hexagon H a)b \ (> k U<fc). By enumerating 
tilings of H ay b \ (>& U <*,) in terms of this factorisation formula one obtains, for free, an 
enumeration formula for horizontally symmetric tilings of H a b \ (> k U <^). 

A similar correlation function involving H+ b \ (t> fe U <^) may also be defined, namely 


u H +(k]t,) 


lim 

a—>oo 


MjHjj, \(> t U < t )) 

M (Hit) 


which gives the interaction between a pair of triangular holes of side length two that lie 
on the fixed weighted boundary of a sea of unit rhombi. According to [12], the number 
of rhombus tilings of Hf h is equal to the number of tilings of V (l)b and so this correlation 
function becomes 


un+ik-f) 


lim 

a—>oo 


M(H+ b \ (>„ U <»)) 

M( V*,i) 
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Theorem 12.21 above is obtained by directly analysing ov(/c;£) as k becomes large, 
whereas Theorem 12.11 follows from considering the product of u H +[k ;£) and u H -[k]ff) 
asymptotically. In order to do so, it is necessary to first find exact enumeration for¬ 
mulas for M(V a . h \ > fc ), M(H~ b \ (> k U < fc )), and M(H+ b \ (> k U < k )). Sections H U 
and [5] are concerned with stating and proving these exact formulas, while the proofs of 
Theorems 12.1112.21 and 12.31 may be found in Section |6j 


3. Some Exact Enumeration Formulas 

The enumeration of rhombus tilings (in various forms and guises) has a fairly long 
history beginning with MacMahon [25] in the early 20th century. Although at the time 
MacMahon was interested in counting the number of plane partitions whose Young 
tableaux fit inside an a x b x c box, there is a straightforward bijection between the 
three dimensional representation of these objects (as unit cubes stacked in the corner 
of an a x b x c box) and two dimensional unit rhombus tilings of H a b<c , the hexagon 
with side lengths a, 6, c, a, b, c (going clockwise from the southwest edge). Under this 
same bijection symmetric plane partitions whose Young tableaux fit inside an a x a x b 
box correspond to vertically symmetric rhombus tilings of H a b- Similarly, transpose- 
complementary plane partitions whose Young tableaux fit inside an a x a x 2b box 
correspond to horizontally symmetric tilings of H a 2 b- 

In |25l p. 270] MacMahon conjectured a formula for the weighted enumeration of 
symmetric plane partitions of height at most b whose Young Tableaux fit inside an 
a x a x b box. Under some specialisation, MacMahon’s conjectured formula counts 
vertically symmetric rhombus tilings of H ab . Andrews managed to prove MacMahon’s 
conjecture in full in [T], and hence established the (ordinary) enumeration of symmetric 
plane partitions. As far as ordinary enumeration goes, however, Andrews was not 
the first to prove the corresponding result. It was Gordon who, around 1970 (though 
published much later in Id). proved the so-called Bender-Knuth conjecture [2], which 
for a special case gives the ordinary enumeration of symmetric plane partitions. Since 
then further refinements and alternative proofs have been found by Macdonald [24l 
pp. 83-85], Proctor [27, Prop. 7.3], Fischer [H] and Krattenthaler J23[ Theorem 13]. 
This theorem states (in an equivalent form) that the number of vertically symmetric 
rhombus tilings of H a b (that is, the number of rhombus tilings of the region V a _f) is 


ST(a,b ) 


n 


2i + b — 1 
2i-l 


n 


i+j+b—1 
i+j-1 


A similar result holds for rhombus tilings of V a ^ \ > k . 

Theorem 3.1. Suppose k is a non-negative integer. 

(i) The number of vertically symmetric tilings of H n ^ m \ (>* U < k ) is 


^ ^ ^n,fc('S) • D r i tk (s') 
s=l 


x ST (n, 2m), 
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where 

(-1 ) s+1 (-k + n + l)!(n + s - l)!(n + 2 s - 1 )! (f + f + s - 2 )! 
= (s-l)!(f-|)!(| + f-l)!(2n + 2s-l)!(—l + f + s)! : 

(-l) s+1 (2s - 2)!(n - k)\(n + s - 1)! (f + f + s - 2)! 
n ’ fc(S) “ (s-l)!(f-|)!(| + f-l)!(n + 2s-2)!(—f + f + s)! 

2 (—l) s+1 (2s - 2)!(n - k)\{n + s)! (| + f + s - 2)! 

+ (» - 2)! (f - |)! (| + |)!(» + 2s - 2)! (-§ + § + s)i ' 

(ii) The number of vertically symmetric tilmgs \ (>& U<fc) is 


ST(n, 2 m — 1) x 


m— 1 


I TP* 7 \ ^n( S )^n,k( 7Tl y 1 D* 7 \ 

T. ( KAA - D . (m) ) ■ KA°) 


S=1 


where 


KAi) = 


2(—1 Y +1 (j + n — 1)!(2 j + n — l)!(n — k)\ (|(2j + k + n — 5))! 


+ 


K(i) = 


K,S) = 


(.j - l)!(2j + 2 n- 1)! (| (—k + n - 1))! (| (k + n- 3))! (|(2j - k + n + 1))! 

4(—1 y +1 (j + n)!(2j + n — l)\(n — k)\ (|(2 j + k + n -* 5))! 

(.j - 2)!(2j + 2 n — 1)! (|(—k + n - 1))! (|(/c + n - 1))! (|(2j — k + n+ 1))! 
2 n (2z — 2)!(z + n — 1)! 

(z — l)!(2z + n — 2)! ’ 

(-l)* +1 (2z - 2)!(z + n — l)!(n — A;)! (|(2z + k + n - 3))! 

(i - 1)!(2 i + n — 2)! {\{-k + n - 1))! (\{k + n — 1))! (|(2 i — k + n — 1))! 


The following formula, due to Proctor [27], counts the number of transpose-complementary 
plane partitions that fit inside an a x a x 2b box: 


TC(a, 2b) 


fa + b — l\ I r I r 2b + % + j + 1 

V a- 1 ) Al-U i + j + 1 

X 7 Z— 1 J 


Since transpose complementary plane partitions correspond to horizontally symmetric 
tilings of H a: b, the above formula counts precisely the number of rhombus tilings of the 
region H~ 2b . Again, a similar result holds for rhombus tilings of H~ 2b \ (>*, U <+). 


Theorem 3.2. The number of horizontally symmetric tilings of H n ^m\ (t>fcU<+) (equiv¬ 
alently, the number of tilings of the region H~_ x 2m+1 \ (> fe U <*,) J is given by 


where 


B'„A S ) 

KM 


£ka s )-km 

, s=1 


• TC{n , 2 m), 


(—l) s+1 (n — k)\(n + s — 2)!(n + 2s — 1)! (|(k + n + 2s — 4))! 

2(s - 1)!^! (| (k + n - 2))!(2 n + 2s - 3)! (-| + f + s)! ’ 
(—l) s+ 1 (2s)!(n — k)\(n + s — 1)! (|(/c + zz + 2s — 4))! 

2 (s!)^! (|(fc + n - 2 ))!(n + 2 s - 2 )! (-| + f + s)! ' 
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Remark 3.1. Observe that if b is odd there exist no horizontally symmetric tilings of 
the region H ab . 


The following enumerative result counts the number of tilings of the weighted region 
Hf b \ (>fc U <fc) defined in the previous section. 


Theorem 3.3. The number of rhombus tilings of Hf 2m \ (> fc U < k ), the upper half of 
the region H n2m \ (>k U <*,) containing lozenges of weight 2 lying along the horizontal 
symmetry axis is 

m 


^ ' B nk (^) ' E nk (^s) 


ST(n, 2m), 


where B n ^(s) is defined as in Theorem \3.1\ and 


(— l) s+1 (2s - 2 )\(-k + n + l)!(n + s - 1 )! (f + f + s - 2 )! 
( S -l)!(f-|)!(| + f-l)!(n + 2 S -2)!(-f + | + S )! ' 


Similarly, the number of tilings of Hf +12rn _ 1 \(> k U<fc) is equal to 2-M(Hf 2m \(> k U< k )). 


Remark 3.2. By re-writing the above result in terms of Theorem 12.21 and combining it 
with Theorem 13. 2 1 via Proposition 12.41 one may construct a factorisation for H ab \ (>*. U 
<fc) that is in a sense similar to that of Ciucu and Krattenthaler [12]. Indeed, it was the 
observation of this factorisation that lead to the results presented in the current work. 


MacMahon proved in [25] that the total number of plane partitions of an a x b x c 
box (equivalently, the number of rhombus tilings of the hexagon H ab( f) is given by the 
formula 


T(a,b, c ) = nun 


i=1j=1k=1 


i + j + k — 1 
i + j + k — 2 


Combining Theorem 13.21 with Theorem 13.31 via Proposition 12.41 gives a similar formula 
for enumerating tilings of H ab \ (> k U <*,). 


Theorem 3.4. The total number of tilings of the region H n ^ m \ (>& U <^) is 




3=1 


X 


■ E n>k (t) 


t= i 


x T{n, 2m, n), 


whilst the number of tilings of H n 2m -i \ (>k U < k ) is 


y ^ B n —l ,fc (^) ' , k (,s) 


S= 1 


x 


771—1 


E B Li ■ £>;.+«(«) 


t= 1 


x T{n, 2m — 1, n). 


Remark 3.3. If k — n, then the hrst formula in the above theorem gives the number 
of rhombus tilings of the hexagon H n 2m containing two unit triangular dents in each 
vertical side, positioned symmetrically with respect to the horizontal symmetry axis. 
This is a special case of the more general class of rhombus tilings considered by Ciucu 
and Fischer in [11 j. 


The following two sections are dedicated to proving Theorem 13.11 Theorem 13.21 and 
Theorem 13.31 by translating sets of rhombus tilings into families of non-intersecting 
paths (Section S]). The enumeration of these families of paths can be expressed as the 
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determinant (or Pfaffian) of certain matrices. These determinants and Pfaffians are 
then evaluated in Section 0 

4. Non-intersecting Lattice Paths 

The aim of this section is to express the number of rhombus tilings of V a j, \ t> fc , 
Hf b \ and H~ b \ (> fc U<fc) as determinants (or Pfaffians) of particular matrices. 

Consider first the region V a ,b\>k described in Section [2j Tilings of this region may be 
represented by unique families of lattice paths across dimers (see Figure 01 left), which 
may in turn be translated into unique families of non-intersecting lattice paths consisting 
of north and east unit steps on the integer lattice with origin O' (see Figure |4] right). 
Counting rhombus tilings of V a ,b \ >k is therefore equivalent to counting the number of 
families (Pi, P 2 ,..., Pf) of non-intersecting lattice paths between a set of starting points 
A = {(— s,s) : 1 < s < b} and a set of end points / = {(x,y) E Z 2 : x + y = 
a} U {t - , t + }, with the requirement that t + = ((a — b — k)/2 — 1, (a + b — k)/2 + 1) and 
t~ — ((a — b — k)/ 2 , (a + b — k)/ 2 ) must be included as end points. Figure Q] shows the 
unique family of lattice paths corresponding to the tiling found in Figure [2J 

By the well-known theorem of Stembridge [HI], such families of non-intersecting paths 
may be expressed in terms of a Pfaffian , defined below. 

Definition 4.1. The Pfaffian of a 2n x 2n skew-symmetric matrix A is 

Pf ( A ) = s § n ( 7r ) II ( A )«> 

7rGA42n i<j matched in n 

where M. 2 n denotes the set of all perfect matchings on the set of vertices {1, 2 ,..., 2 n} 
and sgn(ir ) = (— 1) cr '( 7r ), where cr(n) is the number of crossings of it. 

Remark 4.1. It is a well known fact that 

Pf(A) 2 = det(A). 

Theorem 4.2 (Stembridge). Let V = {rq, ...,ry}, U = {mi, ...,u p } and E = {ei, ...,e g } 
be finite sets of lattice points in the integer lattice Z 2 and assume that p + r is even. 
Denote by PA [A —>■ B) the total number of families of non-intersecting paths from the 
set A to the set B, whilst P(a —> b) denotes the number of paths starting at point a and 
ending at point b. If E is totally ordered in such a way that U is D-compatibl ^ with 
V © E (disjoint union), then the total number of non-intersecting paths starting at U 
and ending at v © E is 

3»(U^,vS>E)= Pf(^ t qV 

where Q = (Qi,j)i<i,j< P is given by 

Qi,j = 'y ] (P(uj —> e s )P(uj —y ef) — P(ui —> et)P(uj -A e s )), 

l<s<t<q 

and Y = (T ij -) ie{lv .. )P}Je{1> ... ir} is given by 

Y^j P [Ui ^ Vj ). 

2 This is a technical condition on the order in which the start points are connected to the end points 
in each family of paths. 
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FIGURE 4. Left: the tiling from Figure [2] as lattice paths across dimers. 
Right: the translation of this tiling to a family of non-intersecting paths 
in Z 2 . 


Remark 4.2. In case p + r is odd, a phantom vertex u p+ i = e q+l may be adjoined with 
the property that 


P(u s —> e q+ i) 


1 , s = p +1 , 
0 , otherwise. 


Then it may be shown that &{U U u p+ \ —* (v © E) U e q+ \) is given by 


Pf 




where Y = (Yi,j)i<i< P +i,i<j< r is as above with the added condition that Y p+ ij = 0 for 
all j, and Q = (Qij)i<ij<p+i is given by 

Qi,j = ^2 ( P{ui ->• e s )P(uj ->• e t ) - P(ui ->• e t )P(uj ->■ e s )). 

l<s<t<g~hl 


Remark 4.3. Throughout this article sums in which the starting index is not necessar¬ 
ily smaller than the ending index are interpreted according to the following standard 
convention: 


n —1 


E"=m Exp(r) 


Exp(r) = 0 


n > m 
n = m 


(4.1) 


- YJr=n Exp(r) n < m. 


Theorem 14.21 may be used to count families of non-intersecting paths that correspond 
to vertically symmetric tilings of holey hexagons. 









































THREE INTERACTIONS OF HOLES IN TWO DIMENSIONAL DIMER SYSTEMS 


13 


Proposition 4.3. The number of vertically symmetric tilings of the hexagon H n2m \ 
(>k U<fc) is the Pfaffian of the skew-symmetric matrix F = (/ij)i<ij< 2 m +2 defined by 


( ( 2n \ 
2^ir=i— - 1 1 ' 1 1 


fi,j ~ < 


■j +1 \n-\-rJ 5 
n—k \ 

(n— k)/2—m-\-i) 5 
/ n—k \ 

\(n— k)/2—m— 1+z/ 5 

lo, 


1 < i < j < 2m, 
i G {1, 2m}, j = 2m + 1, 
i G {1,..., 2m}, j = 2m + 2, 
2m + 1 < i < j < 2m + 2, 


whilst the number of vertically symmetric tilings of H n 2 rn -\ \ (>*, U <*.) is the Pfaffian 
of the skew-symmetric matrix F* = {f*j)\<i. 3 < 2 m +2 defined by 


' ( 2n \ 

2_^ r =i—j+l \n+r) ’ 


f*. = 


n—k 

(n-k-\-l)/2—m-\-i 


n—k 


(n—k—l)/2—m+i 


o, 


1 < i < j < 2 m — 1, 
i G (1,..., 2m — 1}, j = 2m 
i G (1,..., 2m — 1}, j = 2m + 1, 
i G (1,..., 2m — 1}, j = 2m + 2, 
2m + 1 < i < j < 2m + 2. 


Proof. Consider first the region H n2m . \ (>*, U <fc). Let U = A, V = {t“,t + } and 
E = I \ v, and F be the resulting matrix in Theorem 14.21 Then by the simple fact 
that P((a, 6) —* ( c,d )) = ( c_ “^~ 6 ), it should be clear that the proposition holds for 
j = 2m + 1 and j = 2m + 2. For 1 < i < j < 2m, 



n+2m 2ra+n 


= E E 

t=l S=1 



n 

s + t-j 



n 

s + t 


n+2m 


2 n 


^ VV^ + t + f-j 


2n 

n + t + j — i 



where the Chu-Vandermonde convolution has been used in the second line. 

Next consider the region H n2 m-i \ (>khi<k)- Let the point A m+ i = J n+m+1 = (—1, 0) 
be a phantom vertex as described in Remark l4~2l Let U = AUd m+ i, V = {t~,t + } and 
E = (I\v) U/ n+m+ i in Theorem 14 .2 1 and let F* be the resulting matrix. Then F* agrees 
with the matrix F (the determinant of which counts tilings of V n)2m \ >fc) everywhere 
except in the ( 2 m)-th row and column, where 

n+m 

film = P ( Ui ^ es "> 

3=1 

= 2 n , 
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FIGURE 5. Left: A tiling of Hq 6 \ (> 4 U< 4 ). Right: The corresponding 
set of lattice paths. Lines and dashed lines marked in red highlight forced 
dimers on the vertical sides. 


and 

n+m 

fLj = - J2 P ( u i e ^> 

s =1 

= — 2 n , 

for 1 < i, j < 2 m — 1. This concludes the proof. □ 

Consider now the region H~ 2m \(> k U< k ). Observe that in any tiling of H~ 2m \(\> k \J< k ) 
the left and right most vertical sides comprise a set of forced rhombi (Figure O left, 
illustrates such a forcing). Hence it follows that 

l,2m+l \ O ;)) = ^(^n,2m \ O <k)), 

so it suffices to enumerate tilings of the region H~ 2m \ (> fc U < k ). 

Consider then the region H n 2m \ (> k U < k ) for some integer m. Once again, by 
translating sets of tilings to sets of lattice paths across dimers, sets of tilings of H~ 2m \ 
(>fc U <*;) correspond to families of non-intersecting paths beginning at a set of start 
points A = US' As, where A s = (s, 1 — s) for s G {1,..., m} and A m+1 = ((n + k)/2 + 
1, (n + k)/ 2 ), and ending at I = (J ^ 1 I t , where I t = (n +1 , n + 1 — t) for t e {1,..., m} 
and I m +i = {(n — k)/2 + 1 , (n — k)/2), with the condition that no path intersects the 
line y = x (see Figure [5]). With this labelling the number of tilings of H~ 2m \ (\> k U <^) 
is equal to the number of families of non-intersecting paths {P\ , P 2j • • •, P m + 1) where 
Pi denotes a path from Ai to / CT p) and a = (l,m + 1 ) is the permutation on m + 1 
elements. Such families of paths may be counted directly using the well-known theorem 
of Lindstrom-Gessel-Viennot [16], recalled below. 

Let G = (V, E) be a locally directed acyclic graph, that is, every vertex has finite 
degree and G contains no directed cycles. Consider base vertices B = { 6 1; b 2 ,..., b n } 
and destination vertices D = {di, d 2 , ■ ■ ■, d m } and assign to each directed edge e a 
weight, oj e , belonging to some commutative ring. For each directed path P : b —> d, let 
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oj(P) be the product of the weights of the edges of the path. For any two vertices b and 
d. let e(b, d ) = Yp-b^-d 00 ^) denote the sum over all paths from b to d. If we assign the 
weight 1 to each edge then e(b, d) counts the number of paths from b to d. Let V denote 
the n-tuple of paths (Pi, P 2 ,..., P n ) from B to D that satisfy the following properties: 

• There exists a permutation a of {1, 2,..., n} such that for every i, Pi is a path 
running from bj to d a ^. 

• For any i ^ j, the paths P % and Pj have no common vertices. 

Given such an n-tuple, denote by cr(fP) the permutation from the first condition above. 

Consider now the matrix N with (i, j)-entry equal to e(bi, dj). The Lindstrom-Gessel- 
Viennot Theorem then states that the determinant of N is the signed sum over all 
n-tuples V = (Pi, P 2 ,..., P n ) of non-intersecting paths from B to D, 


det(iV) = sgn(a(P))[Ju;(P i ). 

ViB-tD i= 1 

That is, the determinant of N counts the weights of all sets of non-intersecting paths 
from B to D, each affected by the sign of the corresponding permutation a given by 
Pj taking b, to d a (p. Indeed if the weights of each edge are 1 and the only permutation 
possible is the identity then det(IV) counts exactly &(B D ) (that is, the number of 
non-intersecting paths between B and D). 

Proposition 4.4. The number of horizontally symmetric tilings of iL n>2m \ (>& U <*,) is 
(—det(Cr)), where G = ( gi,j)i<i,j<m+i is the matrix given by 


9i,j 


/ 2 n \ / 2 n \ 


n—k \ / n—k \ 

-i) ~ \(n—k)/2—i) ’ 


(n—k)/ 2+1 


/ n—k \ / n—k \ 

\(n— fc)/2+l— j) \(n— k)/2— j) 5 

o, 


1 < f, j < m, 
j — m + 1, 1 < i < m, 
i — m + 1, 1 < j < m, 
otherwise. 


Proof. In the result of Lindstrom, Gessel and Viennot replace B with the set A and 
replace D with the set / (defined previously) and let G be the matrix whose (i, j)-entry 
is the number of paths from Ai to Ij. The only permutation that gives rise to (m + 1)- 
tuples of non-intersecting paths, V m+ i, is the permutation a = ( 1 ,m + 1 ) on m + 1 
objects. By Lindstrom-Gessel-Viennot it follows that det(G) = — &(A —$■ /), where 
PA[X —> Y ) is as in Theorem 14.21 

All that remains to be seen is that the entries of G do indeed count the number of 
paths between Ai and I r Care must be taken to ensure that only those paths that do 
not intersect the line y — x are counted. 

Given a starting point a and an ending point b, let V denote the set of all paths from 
a to b (including those that cross the line y = x). Consider a second set of paths V 
obtained by taking all paths in V that intersect y = x and reflecting portions between 
touching points and the last segment of each path in the line y — x, thereby obtaining 
all paths from a to b' (where b' is the reflection of b in y — x). Then the difference 
\V\ — \P'\ gives the total number of paths from a to b that never intersect the line y — x. 

Suppose 1 < z, j < m. Then the number of paths starting at Ai and ending at Ij is 
given by 

P(Ai -A Ij) = P((i, 1 - i) ->■ (n + j, n + 1 - j)). 





16 


T. GILMORE 


Letting /• denote the reflection of Ij in the line y — x, the number of paths from A, to 
Ij that do not intersect the line y = x is then 


P(Ai -A Ij) - P(Ai -A /') 



n 


2 n 

j -i + 1 


A similar argument holds for the other entries of G. This concludes proof. 


□ 


Precisely the same method may be used to enumerate tilings of the region containing 
weighted lozenges, Hf 2m \ (>*, U <*,), with the exception that any path from A to I 
that intersects the line y — x at f-rnany points has a total weight of 2 t . It is not hard 
to convince oneself that \V\ + \V'\ counts all such paths. Such an argument appears 
in ra- Proceeding in exactly the same way as Proposition 14.41 gives the following: 


Proposition 4.5. The number of tilings of the weighted region H+ 2 m \(>fcU<ifc) is given 
by (—det(G + )) ; where G + = (g£j)i<i,j<m+i is the matrix defined by 

1 ) + L+i'-j) > 1 - i g) < m, 

n + = (mvw) ’ j = m + 1,1 < f < m, 

J 1 tn-lGA’ ! = m+l,l<i<m, 




0, 


otherwise. 


Remark 4.4. Consider the lattice path representation of Hf 2m \(> k U<ik) defined earlier. 
A lattice path representation of Hf l+1 2m _i \ (>fc U <^) may be obtained from that of 
Pn 2 m \ (Pfc U <fc) by extending each start point A s to the left by 1 and extending each 
end point I s vertically by 1 for s G {1,..., m}. Any tiling of this region will contain 
forced lozenges along its vertical sides and so 

M (Hn+i,2m-i \ (>* U < fc )) = 2 • M(F+ 2m \ (> fc U <*)), 
where the extra factor of 2 is due to the starting point A\ lying on the line y — x. 


Theorems 13.1113.21 and 13.31 follow from evaluating the determinants and Pfafflans of 
the matrices defined above. In the following section these determinants and Pfaffians 
are evaluated by considering the Lf/-decomposition of certain matrices. 


5. Evaluation of Determinants 

This section begins with an extension of Gordon’s Lemma m that expresses the 
Pfaffian of the matrix F defined in Theorem 14.31 as the determinant of a much smaller 
matrix. Similar manipulations reduce the matrix F* to the same effect. The (unsigned) 
determinants of these two smaller matrices, as well as those of Proposition 14.41 and 
Proposition 14.51 are then evaluated by finding their unique LI/-decompositions. 

Lemma 5.1. For a positive integer m and a non-negative integer l, let A be the (2 m + 
21) x (2m + 21) skew-symmetric matrix of the form 


A 




for which the following properties hold: 


(i) X is a 2m x 2m matrix such that X = (t./-*) i<*,j<2m and Xj^ = —x t ,j; 
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(%%) Z = (zi t j)i<i t j< 2 i is a matrix satisfying z i} j + z i+ ij + z it j + i = 0 for 1 < i, j < l, and 

%j,i %ij j 

(Hi) Y = (yi,j)i<i< 2 m,i<j< 2 i is a matrix for which 

_ jy 2 m-i,j , 1 < i < m, 1 < j < l, 

\y- 2 m+\-i,j-u 1 < i < 2m,l + 1 < j < 21. 

Then 

Pf(A) = (-l)(a)det(B), 

where B is an (m + /) x (m + /) matrix of the form 



the block matrices of which are defined by 

( X)ij = x i+j -1 + x i+j - 3 H-h x\i-j\+i for 1 < i,j < m, 

i—1 

( Y i)i,j = 5 ~ 2 (ym+i-i+ 2 sj - y m +i- 2 s,j) fori <i <m and 1 <j<l, 

s= 0 

(*2 )i,j = + y m +i-j+ 2 s,i) fori <i< l and 1 <j< m, 

s= 0 

(Z)ij = Zij+i + Zi + i jj+ i for 1 < i, j < l. 

Proof. Beginning with A, construct a new matrix A' = (AY)i< i j< 2 m+ 2 i by simultane¬ 
ously replacing the i-th row of A by the sum 

m—i 

y^(row i + 2s of A), 

s=0 

and the (2m + 1 — i)-th row of A with 

m—i 

y~~](row 2m + 1 — i — 2s of A) 

s=0 

for i = 1, m—1. Perform analogous operations on the columns of the resulting matrix. 
Note that these operations do not change the value of Pf(A) and so Pf(A) = Pf(A'). 
It follows that for 1 < i, j < m, 

m—j m—i 

Am EE %j-i-\-2s—2r 

s =0 r=0 
m—i—j 

= ^ (minjf -f m + 1, m - i + 1} — max{ 0 , t + j})x j+i+2 t- 

t=—m 
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By assumption X is skew-symmetric, so x r = —x r for all 1 < r < 2m. One may also 
verify that 

minjt + m + 1, m — i + 1} 

— max{0, t + j} = min{(— t — i — j) + m + 1, m — i + 1} 

- max{0, (-t-i- j ) + j}, 

whence AY vanishes for 1 < i, j < m. It is not hard to convince oneself that replacing 
i and j with 2m — i + 1 and 2m — j + 1 respectively in the above summation gives the 
(i, j)-entry of A' for m + 1 < i, j < 2m, and so these entries also vanish. 

For 1 < i < m and m + 1 < j < 2m, 

j—m—l m—i 

Ai,j = y ] x j-i-2r-2s 

S =0 7*=0 

j-i 

= y^(rniri{f, m — i + 1} — max{0, t — j + m})xi-j + 2 t- 

t =l 

It is easy to see that 


min{(j — i — t), m — i + 1} 

— max{0, (j — i — t) — j + m} = min{t, m — i + 1} 

— max{0, t — j + m} + 1, 

so again by the skew-symmetry of X, 

Ai t j = x i,ji 

for 1 < i < m and m + 1 < j < 2m, where 

x i,j — x j-i + x j-i-2 + • • • + X\2m-i—j+l\ + l- (5-1) 

Having performed exactly the same operations on both the rows and columns of A 
it follows that A' is also skew-symmetric, thus for m + 1 < i < 2m and 1 < j < m , 
A' = — A' 

For the remaining (i, j)-entries of A' it suffices to consider only those entries for 
which 1 < % < 2m and j G {2m + 1,..., 2m + /}. These columns are affected by the row 
operations alone so for 1 < i < m the (i,j)-entry of A' is 

m—i 

y ^ yi+ 2 s,j— 2 mi 

s=0 

whilst an analogous argument shows that for m + 1 < i < 2m, a' tJ is equal to 

i—in —1 

y ^ Ui—2s, j—2m- 


By the symmetry of Y it should be clear that for j G {2m + l + 1,..., 2m + 21} and 
1 < i < m, 


+i = E 
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whilst for m + 1 < i < 2m and j G {2m + 1 + 1 ,..., 2m + 21} the (i, j)-entry of A' is 

i—m—1 

s =0 

Note that A' has now been completely determined since these row and column opera¬ 
tions leave Z unchanged. 

Construct a new matrix A" = ( A "by perfoming the following operations 
on A'\ 

(i) Add column (2m + l + j) to column (2m + j ) for j G {1,/} (resp. rows); 

(ii) Subtract column (m + j) from column (m+l—j) for j G {!,,..,m} (resp. rows). 

Consider the effect of the first set of row and column operations stated above. For 

1 < i, j < 2 1, 


A" 

‘ rL i-\-2m i j +2 m 


= < 


o, 

%i,j 5 

%i,j Zi,j-\-h 




1 < i,j < l, 

i G {1, •••, /}, j G {/ + 1,..., 21}, 
i G {/ + 1,..., 21}, j G {1, 
l + 1 < i, j < 21, 


by the second property in the statement of the lemma. 

The (i, j)-entry of A" for 1 < i < m and 2m + 1 < j < 2m + l becomes 


-A-ij ^ ^ yi+2s,j—2m T ^ ^ y2m+l—i—2s,j—2mi 
while for rn + !<?'< 2m, 

+ = E yi—2sj—2m ^ ^ ?/2m+l— i-\-2s,j— 2m? 


i—m —1 


i—m —1 


s =0 


s =0 


and crucially (by the third property of the statement of the lemma), A+ = A" m+1 _ lJ 
for 1 < i < m and j G {2m + 1,..., 2m + /}. 

Now consider the effect of the second set of operations applied to A 1 . Again by 
the skew-symmetry of A' this leaves all entries unchanged except those for which i G 
{1,..., m} and j G {2m + 1, ...2m + 21}. For 2m + 1 < j < 2m + l it should be clear 
that this entry vanishes, hence the resulting matrix has the form 


A" 

where for 1 <i,j < m, 


( 

0 

X 

0 

Y ^\ 


- a : 4 

0 

>2 

y-i 


0 

1 2 

0 

Z\ 


~Y} 

—Y t 

1 3 

-z\ 

Z‘2 ) 


i,j %i,j+m 
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is given by (15. ip and 


m—i 

(^1 )i,j = ^ 'j (y2m+l—i-2s,j Vi+2s,j) 
s =0 

1— 1 

(Y 2 )ij ^ ?/m+l— i+2s,j) 

s =0 

2 — 1 

«)« = £ Vm+l—i-\-2s,j 
s =0 

\^1 )zj ”1“ -^i+Zj+Z 

(■^2 )i,j — ^d+ZJ+Z 


for i e {1,.. • ,m}, j e {1,... ,1}, 
for i e { 1 , E { 1 , 

for i e { 1 , E { 1 , 

for 1 < i, j < l, 
for 1 < i, j < l. 


By rearranging rows and columns in exactly the same way A" may be brought into the 
form 


/ 

0 

0 

X 

Y 1 \ 


0 

0 

—Y* 

1 2 

Zi 


-A* 

>2 

0 

y 3 

V 

-Y} 

-z\ 

I 3 

Z 2 J 


(5.2) 


Since the same operations have been performed on both rows and columns this leaves 
the Pfaffian of A" unchanged. By the well-known identity 


Pf 



qJ = (—!)(“) det(P), 


where P is an arbitrary n x n matrix, it follows that 


Pf(A) = Pf(A") 


(-1 


n+l\ 

2 > det 




Reversing the order of the rows 1 tom and multiplying the last l rows and columns by 
—1 gives 


PfG) 


(_l)((" +i )+(?))det 


(—!)«> det (| 


id 





where A", Y\, Y 2 , and Z are exactly those blocks asserted in the lemma. 


□ 


Remark 5.1. Gordon’s Lemma may be recovered from the previous result by setting 

1 = 0 . 


Applying Lemma [5H] directly to the matrix F = (/i,j)i<i,j< 2 m +2 from Proposition 14.31 
results in the following expression for the signed Pfaffian of F, 

Pf(.F) = (-1)0 det(F), 


(5.3) 
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where F = (fi,j)i<i,j<m+i is the matrix with (i, j)-entries given by 


f fi+j-l + fi+j -3 +-b f\i-j\+l, 


1 < i,j <m 


n( .fm+l— i+2s.2m+l fm+i— 2s,2m+i), 1 < i < m, j = m + 1, 


fij — { 


^-n(/,+ m ,-2.s. 2m +l ~b fm+l — j+2s,2m+i) i i — m + 1,1 < j < m, 

[0, otherwise. 

Construct one final matrix F = (/, :j ) L<jj< m +i by performing one last set of row and 
column operations on F, namely for i G {1, ..,m — 1}, subtract the (m — i)-th row from 
row (m + 1 — i), and once again perform analogous operations on the columns. 

For 1 < j < m, 

fm+l,j fm+l,j fm+l,j —1 

fm+j,2m+l "b fm+l— j,2m+l 

n — k \ / n — k 

f n - k )/2 + j) + V(n - k )/2 + 1 - j 
n — k + 1 


(n-k)/2 + j)' 

which agrees with c/^ +1 • (the entries of the matrix whose determinant counts weighted 
tilings of Hf 2m \ (>k U < fc )) for 1 <j< m. 

For 1 < i, j < m, 

fi,j fi.j fi—l,j \fi,j—1 fi—l,j—l) 
fi+j —1 fi—j fi+j—2 ~h fi—j+1 


n 


2 n 

i-j + 1 


2 n 
n + i — j 


which again agrees with gf- for 1 < i,j < rn. 

The remaining entries of F are given by 

fi,m +1 fi,m +1 fi—l,m+l 

= fm+l—i ,2m+l f 

i —1 

"b 2 ^ 'f fm+i+l —2s,2m+l f7Ti — i,~\~2s 

s =1 

n—k \ ( n 

v (n — fc)/2 + 1 — i 


k 


i —1 




s=2—i 


4 (* - 1 ) 


+ T 


2 * - 1 


(n — k)/2 + i 

n — k 

(n — k)/2 — % + 1 


n — k 

[n — k)/2 + s 


(.k — n) i + (n — k)/2 

where the identity ^F™ =0 (^) = (—l) m ( Ar ^ 1 ) has been used in the fourth line. Note 
how none of these operations have changed the value of the determinant, and also that 
fij = gfj for 1 < i < m + 1 and 1 < j < m. 
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The sign in (15.3p may be ignored since these determinants are considered within the 
context of counting families of non-intersecting paths, and so 

M(V n ,2m\>k) = |det(F)|, (5.4) 

where F = (fi,j)i<i,j<m+i is the matrix given by 


A 


f ( 2 n \ | / 2 n \ 

4(i—1) , 2?;—l 

k—n 


+ 


i-\-(n—k )/2 
n—/c+1 \ 

(n— k)/2+j) ’ 


lo, 


n—k 
(n—k )/2 


-i+i/ ’ 


1 < < m, 

j = m + 1,1 < i < m, 
i — m + 1,1 < j < m, 
otherwise. 


Theorem 5 . 2 . The (m + 1 ) x (m + 1 ) matrix F described above has LU-decomposition 

F = L ■ U, 


where L = ( h,j)i<i,j<m+i has the form 

' Aiihj), 

f _ 
iJ ~ 1 , 

0 , 


1 < j < i < m, 
i = m + l,l<j<m 
i — j — m + 1 
otherwise , 


and U = {uij)i<i t j< m +i has the form 

'C n (i,j ), 1 <i<j<m, 

D n k(i), j = m + 1,1 < i < m, 

l,J - ET= 1 B n , k (s) ■ D n , k (s), i — j — m + 1 , 
k 0, otherwise , 


where 

(n)\(i + j - 2)\{2j + n - 1)\ 
n[hJ) ( 2 j - 2 )!(t - j)\(-i + j + n)\(i + j + n — 1 )! ’ 

(-l) J ' +1 (j + n - l)!( 2 j + n - l)!(n - k + l)\(j + (k + n)/2 - 2 )! 
n,k[j ~ (j - l)!( 2 j + 2 n - l)!((n - k)/2)\((k + n)j 2 - l)!(j + (n - k)/2)V 
(n)!(i + j- 2 )!( 2 z + 2 n — 1 )! 

(j — i)\{2i + n — 2)\{i — j + n)\{i + j + n — 1 )! ’ 

(—l) i+ 1 ( 2 f — 2 )!(z + n — l)!(n — k)\(i + (k + n)/2 — 2 )! 
n ’ k ^ ~ (i - l)!( 2 f + n- 2 )!((n - ifc)/ 2 )!((fc + n )/2 - l)!(i + (n - ifc)/ 2 )! 

2 (—l)* +1 (2f - 2 )!(i + n)!(n - Jfc)!(i + (k + n )/2 - 2)! 

”*"(* — 2 )!( 2 i + n — 2 )!((n — k)/2)\((k + n)/2)\(i + (n — k)/2 )\’ 

Proof. Before embarking on the proof proper observe that the (m + l,m + l)-entry of 
(L ■ U ) is 0, whilst for (i, j) e {1, ..., m} 2 , 

A(i)i) ’ C'nifi j) AiC/A) ' Ai(jA)- 
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The proof then amounts to showing that for 1 < i,j < m + 1, 




S=1 


where t = min{i, j}. 

This is equivalent proving that the following identities hold: 

(0 EU Mi, <0 • C„( S , j) = +1 ) + ( n +i-j) for 1 < < < 3 < 


m. 


(ii) EU A *(*> S ) ' D n,k(s) = + i+ln-l)/2 ) {(n-k)/2-i+l) 1 < i < Ul. 

(iii) T,l=i B n,k( s ) ■ c n{s,j) = 2 +j ) for 1 <j< m. 

In order to prove a hypergeometric identity of the form 

n 

U ( n , k) — S(n, k ), (5.5) 

fc=i 

it suffices to find a suitable recurrence relation that is satished by both sides of f!5.5f) . 
Such a recurrence for the left hand side of (15.5ft may easily be found using your favourite 
implementation of the Zeilberger-Gosper algorithm (see, for example [26]). Then veri¬ 
fying that the right hand side also satisfies such a recurrence and checking initial values 
reduces to routine computations. 

Identity (i) above satisfies the following recurrence 

i 

(-(* — j — n))(i +j-n- 1 ) ^(i4 n (i, s) ■ C n (s,j )) 

S = 1 

i +1 

+ {i 2 + i-j 2 +j - 2n 2 - n) ^(Ai(* + M) • C n (s,j)) 

S =1 

z+2 

+ (i - j + n + 2)(i + j + n + l) ^2(A n (i + 2 , s) • C n (s,j)) = 0 , 

S = 1 

whilst the second identity satisfies the recurrence relation 

i 

(i + (k — n)/2 - l)(2i 2 + 2i + (k - n)/2) ^ A n (i, s) ■ D n ^ k (s) 

S =1 

i+1 

+ (2z 2 - 2i + (k - n)/2){i + (n — fc)/2 + 1) ^ A n (i + l,s) ■ D Ujk (s) = 0. 

The third identity satisfies the recurrence 

i 

(.3 + (k- n)/2 - 1) ^(B n , k (s) ■ C n (s,j)) 


S=1 


s=1 


i+i 


+ (j + (n — k)/2 + 1) ^^( B U:k (s ) • C n (s,j + 1)) — 0, 


S=1 


which completes the proof. 


□ 
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Proof of Theorem \3. 1\ part (i). It requires very little work to see that in the previous 
theorem A n (i,i ) = 1 for i e {1, Hence the matrix L contains Is on its diagonal, 

and so the determinant of F is simply the product of the diagonal entries of U, namely 


n 


(n + 2s) n 
(2s - l) n 


m \ 

- • D n>k (t) J . 

t =i / 


The product on the left hand side above is simply a repackaging of ST(n,2m ) and so 
it follows that 


M(V n> 2 m \ >fc) 


^ ^ Bn,k(s) ‘ Dn,k(s) 

s =1 


x ST(n, 2m, n ). 


□ 


Remark 5.2. By letting fc = 0, the above formula gives the number of vertically symmet¬ 
ric tilings of H n ^m with no rhombi protruding across the vertical line segment of length 
2 that intersects the origin of H a t,. It appears that these kinds of tilings, where lozenges 
are forbidden from crossing segments of lattice lines, have not yet been considered in 
the literature. 

Given that fij = gfj for i e {1,..., m + 1} and j G {1,..., m} (and also / m +i,m+i = 
Qm+i m +1 = 0 ) a little further work yields the following lemma. 

Lemma 5.3. The (m + 1) x (m + 1) matrix G + has LU-decomposition 

G + = L■U + , 

where L is defined according to Theorem 1 5. ill and the matrix U + = (ri^)i<jj< m +i is 
given by 

{ Ui t j, 1 < i < j < m, 

En : k(i), 1 < i < m, j = m+ 1, , 

“ ^n,fc(s) • K,fc(s), i = j = m + 1 . 
where E nj k(s) is defined to be 

(—l) s+ 1 ( 2 s — 2 )!(n — k + l)!(n + s - l)!((fc + n)/2 + s - 2 )! 
n ’ k(S) _ (s - l)!((n - k)/2)!((k + n)/2 - l)!(n + 2 s - 2 )!((n - k)/2 + s )\' 

Proof. It should be immediately clear that the (m + l,m + l)-entry of (L ■ U + ) is 0. 
Given the aforementioned agreement between gf ] and fij it suffices to show that 

i An (*’ S) ' = ((n - k)/2 + i) ' 

A straightforward manipulation gives 

s) • E Uj k(s) B n ^(s) ■ C n (s,i ), 

so summing each side over s from 1 to i gives the result. □ 
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Proof of Theorem 1 3.31 Since this /^[/-decomposition is again unique it follows that 

(n + 2 t) n 


det(G+) = 


^ ^ Bn,k^s) ’ -^ra,fc(®) 


s=l 


X 


n <2 1 ~ i), 


The product on the right hand side of (15.61) may be re-written so that 


M« 2m \ (>„ u <*)) = 


^ ' B n ,k ( , s) ' -^n,A»(®) 


s=l 


x ST(n, 2m) 


= -M(//+ +li2m _ 1 \(> fc U< fc )), 


thus concluding the proof. 


(5.6) 


□ 


Consider now the skew-symmetric matrix F* = (f* t j)i<i,j< 2 m +2 defined in Proposi¬ 
tion 14.31 Clearly this matrix does not satisfy the conditions of Lemma 15.11 so instead a 
new matrix F* = (f*i j)i<i,j< 2 m +2 may be constructed by performing the following row 
and column operations on F*\ 

(i) Replace row i with row (* + 2s), and perform analogous operations on the 

columns. 

(ii) Subtract row (2m + 2 ) from row ( 2 m + 1 ), and perform analogous operations on 
the columns. 

In a similar fashion to the proof of Lemma 15.11 for 1 < i,j < m the (i, j)-entry of F* 
disappears, whilst for 1 < i < m and m + 1 < j < 2m — 1 , 


/*« = E /. 


* 

i-\-2s,j * 


s =0 


For 1 < i < m, = (m + 1 - i)2 n and by observing that f* 2m+ i = 




2m+l 


E 


n — k 


\(n — k — l)/2 — mTiT2s 

s=0 


= f*i 


i,2m+2’ 


Hence for 1 < i < m, the entry f* i 2m+1 
operations have been applied, whilst for m + 1 < i < 2m 


r 




n 


k 


vanishes once the second set of row and column 

1 , 

k 


(n — k + l )/2 — m + i) \(n — k — l )/2 — m + i 


n 


2{m 


n 


n 


k + 1 

k + 1 \(n — A; + l)/2 — m + i 


Since F* is skew-symmetric and exactly the same row and column operations have been 
applied to F* in order to construct F*, the matrix F* has now been determined com¬ 
pletely. By reversing the order of the first m rows and columns, and then interchanging 
rows and columns in the correct way, F* may be brought into the form 


0 

-P T 


P 

* 
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where P = (Pi,j)i<ij< m +i is the matrix with entries given by 


P 


hJ 


( 


< 


E i— 1 f* 

s= 0 «/ m+l—i+2s,j+m? 


2 

y l_ 

Z—/s=i 


i-1 


n—k 


s=0 \(n—fc+l)/2—i+2s 
2j / n—/c+1 

fc+l)/2 —j 


). 


/ n—/c+1 \ 

n—/c+1 \(n—/c+l)/2— j) ’ 


,* e G {1,m — 1}, 

j = m,ie {1, •••, m}, 
j — m + 1,1 < i < m, 
i — m + 1, j G {1,..., m}. 


Then clearly |Pf(F*)| = |det(P)|. Proceed by constructing a final matrix F* = 
(fi,j)i<i,j<m +1 from P by subtracting row (m — i) from row (m + l — i) for i G 
{1,... ,m — 1}, and similarly for j G {2,... ,m — 1} subtract column (m — j) from 
column {m + l — j). 

For i G {1,..., m}, f* m = 2 n {% + 1 — i) = 2”, while 


,, G. / n — k \_ + ( n ~ k i 

‘■ m+1 “ \(n — /p + l)/2 — i + 2s/ +f(n-fc + 3)/2-i + 2d 

/ n — k \ 

\(n - k + l)/2 - y ’ 

where the identity y™ =0 (—l) fc (^) = (—l) m ( A m 1 ) has been used in the second line. For 
1 < j < m — 1, 


/, 


2j 


ra+l,j 


n — A; + 1 


2(j - 1) 


n — k + 1 \(n — k + l)/2 — j J n — k + 1 \(n — k + l)/2 — j + 1/ ’ 


n — k + 1 


and also /^ +1>m = P m +i, m - 


Lastly for i G {1,..., m} and j G {1,..., m — 1}, 

fi,j = (Pi,j — Pi,3~ l) — {Pi-1,j — Pi-1,j-l) 


i— 1 


E</» 


— f* ) 

m+l—i+2s,j+m J m+l—i+2s,j+m—l) 


s=0 


i—2 




^ ^ (•/'m+2—i+2sj+m ./m+2—i+2s j+ra—lJ 

s=0 

2 n 


2n 

n + i — j J \n — i — j + 1 


= fi,r 


Thus it follows that 

M{Vn t 2m- 1 \ >fc) = | det(P*)|, 
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where F* = (f*j)i<i,j<m+i is the matrix given by by 


[ fi,j > 
2 n 


/ n—k \ 

) \(n—k+l)/2—i) ’ 


fh = < 


i G G {1,,.., m — 1}, 

j = m,ie {1,. • • ,m}, 
j — m + 1,1 < i < m, 


11_( n ~ k + 1 t _ gozll f n - fc+1 I i m 4-1 7 g 11 m - 11 

n—k+l \(n—k+\)/2—j) n-k+l \(n-k+3)/2-j) ’ 1 7 t |i, 1/, 

— 0 ?T^wn ) , i = m + 1, j = m 

n—k-\-l \{n—k-\-l)/2—mJ ’ ’ ^ 


lo, 


otherwise. 


Theorem 5.4. T/ie matrix F* defined above has LU-decomposition 

F* = L* ■ U* 

where L* = is defined as 

'A n (i,j), 1 <j<i<m, 

Bfik(j), i — m + 1,1 < j < m 

( _ SsT/ D*(s) ■ Bfi k (s)) , i = m + 1, j = m 


l*. — / 
l i,3 ' 


1 

and U* = (u* j)i<i,j< m +i is given by 


i = j = m + 1 
otherwise. 


C n (i,j), 

1 < 

i < j < m — 1, 

D* n (i), 

j = 

m, 1 < i < m, 

E*, k (i), 

J = 

m + 1,1 < i < m 

P n,k( m ), 

i = 

j — m + 1, 

o, 

otherwise, 


with A n (i,j) and C n [i,j ) as before and 

R * _ 2(—l) J+1 (j + n - l)!(2j + n - l)!(n - fc)! (|(2j + fc + n - 5))! 

(j - 1)!(2j + 2n - 1)! (|(-fc + n - 1))! (\(k + n- 3))! (|(2 j - k + n + 1))! 
4(—l)- 7 " 1 " 1 ^’ + n)!(2j + n — l)\(n — k)\ (|(2 j + k + n — 5))! 


+ 


^, fc (0 = 


(j - 2)!(2j + 2n - 1)! {\{-k + n- 1))! (§(fc + n- 1))! (\(2j - k + n+ 1))! 
2 n (2i — 2)!(i + n — 1)! 

(i-l)!(2i + n-2)! ’ 

(—l)* +1 (2i - 2)!(i + n - 1 )!(n - A;)! (|(2i + k + n - 3))! 

(i - l)!(2i + n - 2)\ (|(-fc + n — 1))! (\{k + n — 1))! (|(2 i - k + n- 1))! 


771—1 


= X! 


S=1 


AiM 


-^(s) 


Proof. It should be clear that the (i, j)-entry of L* -U* is /* for (i,j) G {(m+1, m), (■ m+ 
1 , m + 1)} and since f* :) = fij for (i,j) G {1,... ,m} x {1,..., m — 1}, proof of this 
theorem reduces to showing that the following identities hold: 
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(i) EJ =1 A n{h s) ■ D* n (s) = T for i G {1,..., m}; 

(ii) EUi s) ■ E* n k {s) = ( (n _E“iV 2 -i) for i G {1,..., m}; 

GiO R* fs'l.r' (o — 2 t ( n-k+l \ 2(J—1) ( n-k+1 \ 

' ' 2-^ds=\ n,k\ ) n\ n— fc+1 \(n— fc+l)/2— j) n—k-\-l\(n—k-\-3)/2—j) 

!}• 


for j 6 {1,... m — 


Identity (i) may be written as a hypergeometric series (an explanation as to how this 
may be done is postponed until Section [6]) 


Y,A„(i.s)-D' n (s) 

S=1 


4 r 3 


n + 1, (n + 3)/2, 1 — i, i 
(n + l)/2, n + l + i, n — i + 2 



x 


2"(n + l)(n!) 2 
(n-i + l)!(n + i)!' y ’ 

By applying the following summation formula to (15.7ft (see [29], (2.3.4.6); Appendix 
(III.10)]), 


T(a — b + l)T(a — c + 1) 

T(a + l)T(a — b — c + 1) ’ 
the right hand side of (15.71) may be transformed into 

2 n (n + l)(n!) 2 T(n - i + 2)T(i + n + 1) _ 

(n — i + l)!(n + i)\ r(n + l)T(n + 2) 

Proof of the remaining identities follows in much the same way as Theorem 15.21 The 
left hand side of identity (ii) satisfies the following recurrence 

i 2+1 

(2i + k - n - 1) s) ■ E* k (s ) + (2i - k + n + 1) ^ A n {i + 1, s) ■ E* k (s) = 0, 


4-C3 


a, a/2 + 1, b, c 
a/2, a — b + 1, a — c+1’ 


whilst the left hand side of identity (iii) satisfies 


3 

(2 j + k-n- 3) (4 j 2 + 4j + k - n - 1) ^ B* nk (s) ■ C n (s,j)+ 

3=1 

i+i 

(4 J 2 - 4j + k - n - l)(2j - k + n + 3) B* nk {s) ■ C n (s,j + 1) = 0. 

3=1 

After verifying the right hand side of these identities also satisfy the respective recur¬ 
rence relations and checking initial conditions the proof is complete. □ 


Proof of Theorem \3.1\ part (ii). Again it is easy to see that l* t = 1 for i e {1,..., m+1} 
so it follows that 


det(F*) 


m \ 

JJC+s+j • D*{m) ■ P* k (m). 


The product above may be re-written as 


M(V n , 2 m -1 \ >fc) = Pn,k( m ) x ST(n, 2m - 1), 
thus completing the proof of Theorem 13.11 


□ 
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Remark 5.3. Note that the product in Theorem 13.11 is a re-packaging of the formula 
for the number of symmetric plane partitions in a n x 2 m box. The appearance of 
MacMahon’s original formula for vertically symmetric tilings in Theorem 13.11 is not 
surpising if one considers the translation of such plane partitions to non-intersecting 
lattice paths. After a little thought one sees that these may be counted by the deter¬ 
minant of the (m x m) submatrix obtained by taking rows (resp. columns) 1 ,... ,m 
of the matrix F (or indeed F*) defined above. The determinant of this sub matrix is 
simply IT = i Ui,i = nil! u h = ST ( n > 2m )- 

What remains is to find an explicit formula for the determinant of the matrix G from 
Proposition 14.41 The LU- decoiuposition of this matrix follows below. 

Theorem 5.5. The matrix G whose determinant counts the number of horizontally 
symmetric tilings of H n 2 m \ (>* U <&) has LU-decomposition 

G = L-U , 

where L = (h,j)i<i,j<m+i is given by 

r K(hj), 1 <j<i<m, 

B ' n ,k(j), i = m + 1,1 < j < m, 

1, i = j = m + 1 

0, otherwise , 

and U = (ui,j)i<i,j< m+ i is given by 


h.j < 


u. 


Cnfrj), 


hj 


= 


D 


n,k V 


1 < i < j < m, 
j — m + 1,1 < i < m, 


- E”, B 'n,k ( S ) • D n,k( S ), i= 3 = m + 1 , 

0 , otherwise. 


where 


A 'n(hj) = 

KkU) = 
C' n (i,j) = 


n\(i + j - 2)!(2 j + n- l)!(2i - 1) 

(2 j - l)!(i - j)\(j - i + n)\(i + j + n - 1)! ’ 

(-1 ) J -\j +n- 2)\{2j + n - l)!(n - fc)! (j + | + f - 2)! 

2 (J - l)!(2j + 2n - 3)! (f - |)! (| + f - 1)! (j - | + f)! ’ 
(2 j - 1 )n\(i + j - 2)!(2 i + 2 n- 2)! 


(j — i)\(2i + n — 2)\{i — j + n)\(i + j + n — 1)! ’ 

_ (—l)* +1 (2i)!(i + n — l)!(n — k)\ (i+ | + f-2)! 
n ^ l) ~ 2{i\)(2i + n - 2)! (| + f - 1)!^! (i - | + f)! ' 

Proof. It should be immediately obvious that the (m + l,m + l)-entry of (L ■ U) is 0, 
after a little further thought one sees that 

A 'n(iJ) ■ C' n [i,j) = A' n (j,i) ■ C' n (j,i). 

Straightforward manipulations give 

A 'n{p , Q) ■ D n,k (?) = C' n (p, q ) • B' n k (q), 
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hence it suffices to consider only those entries of L g ■ U g for which l<i<j<m + l. 

The approach here is the same as in the proof of Theorem 15.21 that is, a set of 
recurrence relations are given that satisfy each side of the following identities: 


(0 EU <(<■>) -<Z(sj) = LIU - 

(■■) EU4.M ■ DU*) = ( ( „-*Ui-<) - (<»-A-<)- 

Again by using your favourite implementation of the Zcilberger-Gosper algorithm it is 
possible to show that the left hand side of the first identity above satisfies the following 
recurrence 


(j + n 


i)(i + j - n - 1)^ A'rXh s ) • C' n (s,j) 

3=1 

i+1 

- 2 (i 2 + i — j 2 — j — n 2 - n) Aji + 1, s) ■ Cjs,j) 


i +2 

- (i + j + n + l)(i- j + n + 2)^2A' n (i + 2,s) ■ C' n (s,j) 

S =1 


whilst the second identity satisfies 


0, 


(2* + l)(2i + k - n - 2) A'Ji, s ) • D' n k {s) 

S =1 


+ (2« — l)(2i + n 


i+l 

^ + 2 ) + I; s ) ' D'n,k ( s ) 

S=1 


0. 


Verifying that the right hand side of (i) and (ii) satisfy the corresponding recurrence 
relations and checking intial conditions is simply a routine computation. This completes 
the proof. □ 


Proof of Theroem \3. 6 A Once again the LU -decomposition above is unique, it therefore 
follows that 


det(G) 


-£ B U'0-- D h« 

3=1 


TT i Tl + 2 1 — l) n 


The product on the right hand side may be re-written as Proctor’s formula [30], giving 


U< t )) 




s =1 


x TC(n, 2m) 


M(H -_, 

,2m+l \ (^fc U <lfc))- 


□ 


Remark 5.4. Combining Theorem 13.21 with Theorem 13. 3l bv way of Proposition [2741 yields 
Theorem 13.41 

This concludes the proofs of Theorem 13.11 Theorem 13.21 and Theorem 13.41 The 
asymptotics of the correlation functions defined in Section [2] may now be derived using 
these exact enumeration formulas. 
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6. Asymptotic Analysis 


This final section is devoted to the proofs of Theorems 12.1112.21 and 12.31 The exact 


enumeration formulas from Section [3] are inserted into the corresponding correlation 


functions defined in Section [2] and asymptotic expressions for the different interactions 
are derived. For the sake of simplicity the formulas that enumerate tilings of V a ^ \ >&, 
H+ b \ (> fc U <fc) and H~ b \ (> k U < k ) for even b are used in the corresponding correlation 
functions, though of course the results obtained here agree with the case for odd b. 

Recall the correlation function ujv(k-,£) defined in Section [2l 


where m ~ C.n/2 and M(R ) denotes the number of tilings of the region R. The enumer¬ 
ation formula for M ( 14 , 2 m\>fc) given in Section [3] lends itself readily to asymptotic anal¬ 
ysis of u v {k\ f) since it is of the form M(K, 2 mW) = M( 14,2m) • Y^=i B n,k{ s ) ' Ai,fc(s)- 
It follows that 



Similar arguments hold for uj H - (k]£) and w H + (&;£). The product of these two functions 
then gives U!n(k', 4 according to the Matchings Factorization Theorem [7]. 

In order to extract asymptotic expressions for these correlation functions the finite 
sums from Theorems 13.11 13.21 and 13.31 are first re-written as limits of hypergeometric 
series. By applying some transformation formulas and sending first n and m to infinity, 
one obtains these correlation functions in terms of k (which parametrises the distance 
between the holes). Asymptotic expressions for the different interactions are then de¬ 
rived by letting k become large. Before embarking on the proofs proper, some standard 
definitions and necessary transformation formulas concerning hypergeometric series are 
recalled below. 

Consider the p F q hypergeometric series 



where (ct)p denotes the Pochhammer symbol, that is, 


\ I {p)p — {cx)(a + 1) • • • (a + f3 — 1), 
a )P = 1 


/3 4 0 , 

(3 = 0. 


(«)/3 


A hypergeometric series is considered well-poised if ai + 1 = 4 + <A = • • • = b q + a p . A 
very well-poised hypergeometric series is a well-poised hypergeometric series for which 
o, 2 = ai/2 + 1. 

It is often convenient to re-write a finite sum with summands that consist of hyperge¬ 
ometric terms as a limit hypergeometric series by introducing a term that “artificially” 
terminates the series after a particular point. Consider, for example, the following finite 
sum, 



( 6 . 1 ) 
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By multiplying the summand by (1 — t) s /( 1 — t + e) s it is possible to re-write this 
finite sum as a limit of a hypergeometric series: 


E (0l)s(02)s ' ‘ ' (flp)s Z S 
(bi) s (b 2 )s---(b q ) s s\ 


y ( ( 0 l)s( 02 )s ' ' ' (Op)s(l — t) s 

™ ^ ( 6 i)-( 62 ).-”( 6 9 )-(i-t + e) s ^r; 


lim 

e —>0 


(p+lFq+1 


Oi, 02 ; • • •, Op, 1 — t + e, 

6 1; b 2 , • • • , dq, 1 — t 



Suppose in addition that the terms a±,... ,a p and bi, ... ,b q satisfy the conditions for 
a very well-poised hypergeometric series. Then by multiplying the summand in (16.Tj) 
by a factor of 

(1 - t) s (Oi + t) s 
(a i + t + e) s (l — t + e) s 

the hnite sum (16.If may be expressed as a limit of a very well-poised p+2 F q+2 hyperge¬ 
ometric series 


(ai + c) s (a 2 + e/2 ) s ■ ■ ■ (Op + e/2) s (l — t) s (ai + t) s z s 
e ~l (bi + e/2) s (b 2 + e/2) s ■ ■ ■ (b q + e/2) s (ai + t + e) s (l — t + e) s s! 


lim 


E 


v s=0 

= lim 

e-J-0 


p+2-fg+2 


ai + e, 02 + e/ 2 , 


Op + e/ 2 , 1 — t, ai + f 


bi + e/ 2 , b 2 + e/ 2 , ..., b q + e/ 2 , ai + f + e, 1 — £ + 


There are a wealth of transformation and summation formulas that may be applied 
to hypergeometric series. Throughout this section frequent use shall be made of the 
following transformation formula for a 7 F 6 series (see [29, (2.4.1.1), reversed]) 


7-^6 


a, a /2 + 1 , 6 , c, d , e, —n 
a/ 2 , a — 6 + 1 , a — c+ 1 , a — d + 1 , a —e + 1 , a + n + 1 




(o T l) n (a — d — e T 1) 


(cl — d T l) n (a — e T 1) 


n t? 

- 4-^3 


a — 6 — c + 1 , d, e, —n 
a — 6 + 1 , a —c+ 1 , —a + d + e —n’ 


( 6 , 2 ) 


Another transformation formula which will be useful is the following (see [29, (1.8.10)]) 


2^1 


—n. a 




(1 - 

( c )n 


2 Fi 


—n, c — a 
—a — n + 




-1 


(6.3) 


Proof of Theorem \2.2i Consider the region 14 , 2 m \ »&• According to Section [2] the cor¬ 
relation function of this triangular hole with the free boundary to the right of V n - 2m \>*. 

is 

Bn,k(s') ’ Dn,k(s') I , 


Uv(k;£) = hm 
n—>oo 


m 

E 

V s = 1 


(6.4) 
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where m £n /2 and 

(-1 ) s+1 (-k + n + l)!(n + s - l)!(n + 2 s - 1 )! (| + | + s - 2 )! 
n ’ k{S) = (s - 1 )! (f - |)! (| + f - 1 )!( 2 n + 2s- 1 )! (-| + f + s)\ ’ 

(-l) s+ 1 ( 2 s - 2 )!(n - k)\{n + s - 1 )! (§ + f + s - 2 )! 

“ ( S -l)!(f-|)!(| + f-l)!(n + 2 S -2)!(—| + f + S )! 

2(—l) s+1 (2s - 2)!(n - fc)!(n + s)! (| + § + s - 2)! 

The expression on the right of (j6.4[) may be written as two separate finite sums consist¬ 
ing of purely hypergeometric terms. These finite sums may in turn be expressed as a 
limit of two very well-poised hypergeometric series that have been artificially terminated 
and so (16.4jl may be re-written in the following way: 


k-V+O — hm lim 


n\(n + l)!(n — k)\(n — k + 1)! 
rwooe— >6 \^( 2 n + l)!((n — k)/2)\ 2 ({n — k) /2 + l )! 2 


X 7-^6 


n + l + e, k±n±e, l±i ) n + m + l, l- m 


n+l+e n—fc+4+e n—fc+4+e 2n+3+e 


, 1 — m + e, n + m + l + e 


; 1 


2 ’ 2 ’ 2 ’ 2 

4(n + l)!(n + 3)!((n + k)/2)\{n — k)\(n — k + 1)! 
(2n + 3)!((n - A;)/2)! 2 ((n - k )/2 + 2)! 2 ((n + fc)/2 - 1)! 


X 7 ^6 


o | r n+5+e fc+n+2+e fe+n+2+e 3+e 
° ' e ’ 2 1 2 ‘ 2 '2 
n+3+e n— fc+6+e n—fc+6+e 2n+5+e 


n 


, n + m + 2, 2 — m 


, 2 — m + e, n + m + 2 + 6 




2 ’ 2 ’ 2 ’ 2 
Applying (16.2ft separately to each 7 F 6 series in the above sum gives 


(n + 2 + e) m _i((l + e )/2 — m) m _i 


- 4 -C 3 


(n + (3 + e)/2) m _i(l — m + e) m _i 
for the upper hypergeometric series in (16.5p and 
(n + 4 + e) m _ 2 ((l + e)/2 — m ) m _2 


2 — k, 777 _|_ 77 -)- 1 I— 777 , 

^ ' I 

n—fc+4+e n—fc+4+e 3—e ’ 


+ + (5 + e)/2) m _ 2 (2 — m + e ) m _2 


4-^3 


2 — k, — 

n—fc+6+e n—fc+6+e 5—e 


, m + n + 2 , 2 — m 


(6.5) 


( 6 . 6 ) 


(6.7) 


for the lower. Combining 
tend to zero gives gives 


2 ’ 2 ’ 2 
with the corresponding pre-factor in (16.5ft and letting e 


(n — k)\(n — k + l)!(n + m)\{2m — l)!(n + m — 1 )! 

((n — k)/2)\ 2 ((n — k)/2 + 1 )! 2 (m — l)! 2 ( 2 n + 2 m — 1 )! 

2 — A:, 1/2, m + n + 1, 1 —m 
4 3 [(n- fc + 4)/2, (n-ifc + 4)/2, 3/2 ’ 1 


( 6 . 8 ) 


Applying Stirling’s approximation (n! ~ \[2jvn(n/e) n ') to the pre-factor of (16.81) as n 
tends to infinity (and remembering that m ~ £n/ 2 ) gives 


(etf + a)) 1 /^ 1 -* 
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which implies that (j 6 . 8 j) vanishes in the limit. Hence the correlation function is domi¬ 
nated by the second sum (16.711 . chiefly 


uv(k]Q 


f 2{{n + k)/2)\(n — k)\(n — k + 1)! 

71—> OO V((n - k)/2)\ 2 {{n - k)/ 2 + 2)! 2 ((n + k)/2 - 1)! 

{2m — l)!(m + n + l)(m + n — l)!(m + n)! 

3 (m — 2)!(m — l)!(2m + 2n — 1)! 

2 — k, 3/2, m + n + 2, 2 — m \ 
X4 3 [{n-k + 6)/2, (ra-lb + 6)/2, 5/2’ X J y ’ 


(6.9) 


which follows from letting e tend to zero in (16.711 . 

In order to completely derive the asymptotics of (16.9ft it must be shown the limit and 
sum operations commute in the following expression: 


/\ (2 — k) s {2>/2) s {m + n + 2) s (2 — m) s \ 

™ (((" - k + 6)/2) s ) 2 (5/2) s (s!) ) ■ 

As this series terminates for s > k — 2 , this is in fact a sum with a finite upper bound 
dependent only on k and so the limit and the sum in the above expression may be 
safely interchanged. Doing so, and then applying Stirling’s approximation as n tends 
to infinity, expression (16. 9 p becomes 


uv{k]£) 


(£(£ + 2)) 3 / 2 4 1 ~ fc 
37 t • e 


- 2^1 


2-k, 3/2 
5/2 


-e(e + 2) 


( 6 . 10 ) 


To uncover the asymptotic behaviour of u>v(k]£) as k becomes very large it is first 
convenient to apply transformation (J6j3j) to the hypergeometric series in (16.1011 . thereby 
obtaining 


u v (k]£) 


(£(£ + 2)) 3 / 2 4 1 ~ fc 
37T 


^ 3 (^ + l ) 2fc - 4 F 
V 2{k- 1/2) 2 1 


2 -lb, 1 
3/2 - k 


;(£ + i) 



( 6 . 11 ) 


Consider the hypergeometric series in <ny(A;;£) above, that is, 


(2 - k) s 

A(3/2-fc),(f+l)2'' 


It is easy to convince oneself that for any positive even k, ^ 2 -k) — ^ f° r non_ 

negative s and so 


lim y _ (2 ~ k)s _ 

fc_ 5 , 00 “i) (3/2 ~ k) s (£ + l) 2s 


k—2 


< 


lim 

h —Vnn • ^ 


V (£ + !) 2 


S 


( 6 . 12 ) 


Since £ > 0, the limit on the right hand side of (16.12)1 exists and is finite so the limit 
and sum operations on the left side of (16.12)1 may be interchanged. Hence as k becomes 
very large, 


u v (k-,0 ~ 


a/£(£ + 2 ) 

2nk 


(g + 1) 

2 


2(fc—1) 


which concludes the proof. 


□ 
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Proof of Theorem 13.31 Consider the correlation function 


u H - 


(fc;0 = lim ( 5 Z B 'nA S ) ■ D 'n, 


. S=1 



where again m ~ ^rz/2 and 

, (— l) s+1 (n — k)\(n + s — 2)!(n + 2s — 1)! {^(k + n + 2s — 4))! 

Bn ’ k ^ S ~ 2(s - 1)!+^! (i(fc + n - 2))!(2n + 2s - 3)! (-§ + f + s)! ’ 

, (—l) s+1 (2s)!(n — fc)!(n + s — 1)! (|(/c + n + 2s — 4))! 

Dn ’ k ^ ~ 2(s!)+^! (i(fc + n - 2))!(n + 2s - 2)! (-■§ + f + s)!' 

The above finite sum may be re-written as a limit of a hypergeometric series: 

(l y y i (n - l)!(n + l)!(n - k)! 2 

— lim lim 

n-* * * - *° 

X 7 ^6 


n—»<x>e—>0 ^2(277, — l)!((n — /c)/2)! 2 ((n — fc)/2 + l)! 2 

n + l + e, ^±£, 2±£, n-hm + l, 1-m 


2 ’ 2 ’ 2 
?i+l+e n—fc+4+e n—k+ 4+6 2w+l+e 
2 ’ 2 ’ 2 ’ 2 


’ 2 .. ; 1 

, 1 — m + e, m + n + 1 + e 


Applying transformation formula (16.21) to the hypergeometric series above and letting 
e tend to zero gives 


UH-(k;£) = lim 


(2 m + l)!((n — A;)!) 2 (m + n — 2 )!(m + n)\ 


T, — fc | 2 ( ( 1 
X 4 -F 3 


12(m - l)!m! (^!) {{\{-k + n + 2))!) (2m + 2n - 3)! 

2 — k, 3/2, m + n + 1, 1 


m 


_(n — fc + 4)/2, (n — fc + 4)/2, 5/2 




Consider the limit of this expression as n tends to infinity. By interchanging the sum 
and limit of the hypergeo metric series (as in the previous proof) and applying Stirling’s 
approximation to the pre-factor it may be shown that 


= -3FT- 2Fl 

which is precisely expression (16.10p above. 


'2 — As, 3/2 
5/2 


; -m + 2 ) 


□ 


Proof of Theorem I A 11 By the Matchings Factorization Theorem [7] it suffices to con¬ 
sider the asymptotics of the expression 


^ff+(fc;0 = lim y^B nik (s) ■ E nk (s ) , 

77.—VOO \ L / 


, S = 1 


where of course m ~ fn/ 2, B Ujk (s) is as before and 

(—l) s+ 1 ( 2 s - 2 )!(—A: + n + l)!(n + s - 1 )! (f + f + s - 2 )! 


E n ,kf s) 


(s - 1 )! (f - I)! (| + f - l)!(n + 2 s - 2 )! (-f + f + s)! ' 
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Again this correlation function may be written in as a limit: 

n\(n + l)!(n — k + l )! 2 


i^H+(k:£) = lim lim . . ,... , 

n-t ooe->o \(2n + l)!((n — k)/2)\ 2 ((n — k)/2 + l )! 2 


x?F 6 


71 + l + e, ^±|±i, ±±£, n + m + 1, l-m 

n±+te^ n-k+4+e ^ n-k+ 4+6 ^ 2n+3-H ^ _ m + ^ m + n + 1 + e ’ 


(6.13) 

l , —: 1 - 1 - i — m _i_ c m —i— / >o —i— i —i— i 

2 ’ 2 

Transforming the above hypergeometric series according to (16.21) and letting e tend to 
zero gives 


= lim 


( 2 m — l)!(n — k + 1 )! 2 (m + n)!(m + n — 1 )! 
n->oo \^(m — l)! 2 ((n — fc)/ 2 )! 2 ((n — /c )/2 + l)! 2 ( 2 m + 2 n — 1 )! 

2 — k, 1 / 2 , m + n + 1 , 1 — m 


X4F3 [(n- k + A)/2, (n-k + 4)/2, 3/2 
As n tends to infinity Stirling’s approximation yields 

\/£(£ + 2 ) 

4 fc_ 1 7 r • e 

for the pre-factor while the hypergeometric series reduces to 

"2 - k, 1/2 




2 F 1 


3/2 


;-£(2 + 0 


( 6 . 14 ) 


(6.15) 


(6.16) 


This follows from interchanging the limit and the summand in (I6.14|) , again by a similar 
argument to that found in the proof of Theorem 12.21 Applying transformation 
to (16.161) . the correlation function uj H +[k ;£) becomes 


(£ + i ) 2fc - 4 

u H +{k\(L)= z7 L / ^ 2F 1 


'2 — As, 1 


2{k — 3/2) [ 5/2 -k 

k—2 

£ 


; K +W 2 

{k- 1 


= I ft I 1 \2fc-4 ~ 1 ~ S ) _ 

“2^ + J £z(k-s- 3/2) s+1 (e + l) 2 ^ 


The asymptotic interaction of the holes is given by 

k—2 


lim u> H +(k]£) = lim 

k —^oo k —^oo 


It should be clear that 

k—2 

lim £ 71— 2 ; 

At—V nn ^ ^ I A" — — 


2fc—4 


;(£ + !) 


(fc — 1 - s)c 


E (/c — 1 — s) 6 

(h — .Q — 


^ (* - » - 3 / 2 ) s+1 « + 1 ) 2 * 


( 6 . 17 ) 


k—2 


< 


lim > 

k—>oc ' J 


(e + i ) 5 


k ^°° j^(k-s- 3/2) s+ i(f + l) 2s s=Q 

so since the limit on the right exists and is finite the sum and limit operations in (16.171) 
may be interchanged. Combining (16.171) with (16.151) and letting k become very large 
gives 


u H +(k-,Q 


(£ + l)\ 2fc 2 1 


2nky/S(Z + 2)' 
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Inserting this into Proposition 12.41 along with the result from Theorem 12.31 yields 


u H (k;£) ~ 





This completes the proof. 


□ 
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